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Large N gauge theories with adjoint matter can be numerically studied using lattice techniques. Eguchi-Kawai 
reductions holds for this theory and one can reduce the lattice model to a single site. Hybrid Monte Carlo 
algorithm can be used to simulate this model. One can either perform an exact computation of the "fermionic 
force" or use pseudo fermions as part of the HMC algorithm. The former algorithm is slower than the latter 

C ~ ) but has the advantage that one can work with any real number for the fermion flavor. Some results using both 

£N| algorithms will be presented. 

Ph 

<D 

^ 1. Introduction 

O 

Lattice studies of vector like gauge theories with adjoint fermion matter with the aim of understanding the 
conformal window has recently attracted considerable attention (see [T] and references therein). The gauge 
group is chosen to be SU(N) and the beta function is 

03 

^7* in aM U-4/j 34 - 32/ 3 A 

D 

where a is the lattice spacing, / is the number of Dirac flavors or adjoint fermions and A = g 2 N is the 't 
l— ~~ ' Hooft gauge coupling on the lattice. The first two coefficients in the beta function are renormalization scheme 
I independent. In order to maintain asymptotic freedom, we restrict ourselves to / < 3. The two loop beta 

function has a zero if / = 2 and this has a motivated numerical studies of SU(2) gauge group with two Dirac 
C^i flavors of fermions in the adjoint representation [2H1]- 

A continuum analysis of the theory with adjoint fermions on R 3 x S 1 with periodic boundary conditions for 

fermions in the compact direction shows that the Zn symmetry is not broken in that direction [5] . An analysis 

on S 3 x S 1 also shows a region where the symmetry is not broken [5] . A lattice analysis of the same theory 
Q\ with Wilson fermions indicates that one can reduce the compact direction to a single site on the lattice and still 

maintain the Zn symmetry |S] . This is expected to be the case for / > ^ [3] and for non-zero quark masses 

^- tuna. 

The model on the single site lattice is defined in Sec. [2] We will use one-loop perturbation theory of this 
model to show that it is expected to the correct continuum phase in Sec. [3j A summary of our non-perturbative 
results will be presented in Sec. [I] 

U 

2. The model 

The action on a single site lattice with one flavor of adjoint Dirac overlap fermion is given by 

S = S g + S f . (2) 

The matrices, [i = 1,2,3,4 are elements of the su(N) algebra and conjugate to the four SU(N) gauge 
degrees of freedom, U^; [i = 1, 2, 3,4. The gauge action is 

4 

S g = -bN £ Tr [U^Upi] . (3) 

The lattice gauge coupling constant is b — -?sjj- 1 The overlap fermion action is 

S / = -2/logdettf 0+ . (4) 



1 This coupling is related to the standard lattice coupling, /3, by b = -^j^z- A value of b = 0.35 corresponds to /3 = 2.8 for N = 2 
and p = 6.3 for N = 3. 
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The Hermitian massive overlap Dirac operator is defined by [T31 [H] 



ff (M)=^[(l + /i)75 + (l-M)e(#)], (5) 
where fi G [0,1] is the bare mass and 

4W = ^±^W ± ; p ± = ^, (6) 

factorizes into two disjoint pieces corresponding to the two chiralities. Note that / = 1 in Q is the correct 
result for a single Dirac fermion in the adjoint representation 2 . We can set / to be half integers and simulate 
Majorana fermions but we can also extend / to any real number in Q. 

The function e(H) appearing in ^ is the sign function of the Hermitian Wilson Dirac operator, H. The 
Hermitian Wilson Dirac operator for adjoint fermions is given by 



-3E M 4(^-^) -4 + m + |E^(^ + ^f). 
= (4-^)75-^(^+^1 (7) 



with 



Let $ be a traceless Hermitian matrix and denote one component of an adjoint Dirac fermion on the single site 
lattice. The action of on $ is given by 

V^ = U^Ul Vfe = Ul$U ll . (9) 

One can verify that H is Hermitian in the usual sense: 

Tv^H<P = [Tr&Hf]* = Tr [(H^)^] . (10) 

Therefore = (Jf*)t and in addition it is also true that Tri7$ = if Tr$ = 0. The same is also true for 



3. Weak coupling perturbation theory 

For the single site lattice theory to reproduce the correct infinite volume continuum theory, the center sym- 
metry that takes 

U^^e^U^; k ll = 0,--,N-l (11) 

with /i = 1, 2, 3, 4 independent of each other should not be broken. In the limit of large N, this amounts to 
saying that TrU^ = for all /i which is equivalent to the statement that the eigenvalues of Un are uniformly 
distributed on the unit circle. The single site perturbation theory is given by 

Um = e^D^e-**" ; £>« = e<S tj . (12) 

Keeping 9^ fixed, we expand in powers of a M . The lowest contribution to S g comes from the quadratic term 

in a M [15] and the lowest contribution to Sf comes from setting et M = 0. Each has jv ^~ 1 - ) two by two blocks 
of the form 

oos(0*-0£) sin(^-^)\ 



2 We are assuming that global topology is completely suppressed and one can restrict the theory to the zero topological sector. 
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with 1 < i < j < N. The remaining (N — 1) x (N — 1) matrix is a unit matrix. Therefore, the gauge held 

effectively has (N — 1) zero momentum modes and N(N — 1) non-zero momentum modes of the form e 1 ^^ 8 ^ 
with 1 < i =/= j < N. If 9 1 ^ for a fixed /j, are uniformly distributed on the unit circle and there is no correlation 
between the different /j,, the single site model will correctly reproduce the momentum integral of the infinite 
volume continuum theory. Our aim in one-loop perturbation theory is to study the distribution of 9^. 

The computation of the fermion determinant reduces to a free field calculation at this order and the result is 



S f = -4/^lnA(6' i -^ 

i¥=3 



4(7V-l)/lnA(0) 



(14) 



where 



, is the phase associated with the boundary condition in the /i direction. The eigenvalues, 



±A(p), are two fold degenerate and given by 



A(p) 



\ 



1 + /i 2 1 - /i 2 



2 

2 



in 



(15) 



^sin^-mj +^ iSin 2 ^ 
The complete result from fermions and gauge fields is 



S = Wln 



E 



i 



sin 



— 4/ln A(6»" — 6» J 




-4(iV-l)/lnA(0). 



(16) 



- and the single site model is not in the correct 
and this choice need not be the minimum. 



If / = 0, the minimum of the action occurs when all 9 
continuum phase [15] . If = 0, S 1 ^ — >oo when all 9 1 ^ 

Overlap fermions reproduce the correct continuum behavior by restricting the full Brillouin zone to a physical 
region around zero defined by 



m>2j2 sm 2 y = 2 E sin2 ' 

t 1 M 



(17) 



We cannot set m to be very large since overlap fermions reduce to naive fermions as m — > oo and naive fermions 
on a single site lattice do not reproduce the correct continuum behavior [9]. We cannot set m to be too small 
since we will not cover a substantial region of the Brillouin zone to realize the correct momentum measure. 
Unfortunately ( 16 ) is a function of AN variables and it is not easy to study it analytically and find the 



minimum. One option is to numerically study this function. In order to find the minimum of S, we consider 
the Hamiltonian 



H 



Jek 



■PS. 



(18) 



ft :i 



For large (3, the Boltzmann measure e~ H will be dominated by the minimum of S. We can perform a HMC 
update of the 7r, 9 system to find this minimum [5]. 

A choice for the order parameters associated with the Zfj symmetries is [TS] 



(19) 



If = |, then the Zn symmetry in that direction is not broken. If 9^ are uniformly distributed in a width 
a < 2ir, then 



lim P u = - 
2 



, 2 a 

1 — — sm — 
a 2 



(20) 



The result of the numerical simulation to look for the minimum of S in ( 16 1 with / = 1 is shown in the left 
panel of Fig. [T] It is clear that we will not reproduce the correct continuum behavior if m < 3. It seems that one 
can take m as large as 8. This analysis does not take into account possibilities of correlations in the different 
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directions. Instead of studying correlations in different directions, we compute the correlated action, S c , with 

9„ = "^ji- and compare it to the uncorrelated action, S u , with 6^ = N J where are different permutations 
for different (x. The difference between S u and S c is plotted as a function of to for several values of N in the 
right panel of Fig.[2]for / = f . It shows that the correlated one is below the uncorrelated one for to > 5. Using 
the above two arguments, we conclude that we need to set 3 < to < 5 in order for the single site theory to 
reproduce the correct momentum measure when / = 1. 




Figure 1: The left panel is a plot of P M as a function of m for several different values of N. The right panel shows a 
plot of S u — S c as a function of m for several different values of N. 



4. Non-perturbative results 

We need to verify if the results obtained in one-loop perturbation theory in the previous section remains 
valid in a full non-perturbative computation. The non-perturbative computation is performed using the Hybrid 
Monte Carlo algorithm. Let be traceless Hermitian matrices that are conjugate to the gauge fields, Up. The 
algorithm starts with one choice for Up. Then, we draw according to a Gaussian distribution. 



The equations of motion for Un are 



dU, 



(21) 



Setting ^7=0 results in 



rr dHp 



dS g dS f 



dr 



dr 



(22) 



The derivative of Sf with respect to r is referred to as fermionic force term and is computationally intensive. 
One can derive exact expressions for the single site model [9,. It involves exact diagonalization of H and the 
computational cost grows like A^ 6 . The advantage of computing the fermionic force exactly is that one can 
work with any real value of /. An alternative approach is to use the pseudo-fermion algorithm to compute the 
fermionic force term [16] . This is less computationally intensive but works only for integer values of /. 

The result for massless fermions (we set \l = 0.01 in the numerical simulation) is shown in the left panel of 
Fig. fusing the exact algorithm for the fermionic force. We can see that the theory will reproduce the correct 
continuum limit if to > 3. We set to = 5 and studied the behavior of the model as a function of the quark 
mass. The result is plotted in the right panel of Fig. [2] We see that the model will reproduce correct continuum 
physics even when the fermions are massive. 



4.1. Chiral symmetry 



In order to study whether chiral symmetry is spontaneously broken, we studied the behavior of the low lying 
positive eigenvalues, < Ai < Aa • • ■, of the hermitian overlap Dirac operator. If chiral symmetry is broken, we 
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Figure 2: The left panel is a plot of P M as a function of m in the full numerical simulation with N = 11, b = 7, / = | 
and fi = 0.01. The right panel shows a plot of Pi (this is the one that is broken first) as a function of (j, for several 
different values of N. 



expect a relation of the form 

Zi = N 2 X^{b) (23) 

where the joint distribution of the scaled variables, Zi, are given by some chiral random matrix model j!7| and 
£(&) is the value of the chiral condensate. 

The chiral Random Matrix theory ensemble for a symplectic matrix, C = a^C^, is 

Z = J [dC,]e- t C "] 2 [det H rmt ] S ; H rmt = ( £ (24) 

with C M being a real square matrix. We expect z\ to be the eigenvalues of H rmt with £(&) being the scale that 
relates these eigenvalues to the eigenvalues of H Q . We want to eliminate the scale set by the chiral condensate, 
and we focus on 




The result as a function of / with b = 5, m = 5, fi = 0.01 and N = 11 is shown in the left panel of Fig(3] By 
comparion with chiral random matrix theory, it looks like chiral symmetry is broken for / = and / = 1 but 
not for / > 1. This would be the case if the non-pertubative beta function has a zero for / > 1. 



4.2. Setting the scale 

In order to get some insight into the non-perturbative beta function, we study a lattice scale as a function of 
the lattice coupling. An L x T Wilson loop operator in the \i — v plane is given by 

W(L,T) = U^Ujul L Ul T . (26) 

The eigenvalues, e t6k ; k = 1, • • • , N of this operator are gauge invariant. Let p(9; L, T, b) be the distribution of 
these eigenvalues with 8 6 [— 7r,7r). This distribution undergoes a transition [TS] at N — » oo as the area, LT, 
is changed at a fixed coupling b: the distribution has a gap at it for small areas and it becomes gapless for 
large areas. There is a critical area A c (b) where the gap closes. There is a universal function [TU] describing the 
distribution in terms of the scaled variables derived from A(b) and 9 in the vicinity of A c {b) and 7T. 
Let 



N (y;A,b) = (det(e* + e" = W(L,T))); A = LT. 



(27) 



G 
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Figure 3: The left panel is a plot of r as a function of / in the full numerical simulation with N = 11, b = 5, and 
fj, — 0.01. The right panel shows the lattice scale as a function of the lattice coupling for / = 1. The two loop result for 
/ = and f — 1 are shown for comparison. 



The region close to y = probes close to ir. Let 

N (y; A, b) = Co (A, b, N) + C x (A, b, N)y 2 + C 2 (A, b,N)y 4 + ---. (28) 

It is useful to define 



C ( A,b,N)C 2 (A, b,N 
CKA,b,N) 



Sl(A,b,N) = ^; T^n ■ (29) 



One can show using the universal scaling function that 



r 4 (I) 

ft(A c (b),n, oo) = 48 ^/ = 0.364739936 (30) 



We can define A c (b, N) at a fixed N and b as the area where 

n(A c (b, N),b, N) = 0.364739936, (31) 

and 

lim A c (b,N) = A c {b), (32) 

will be the location of the transition at infinite N. 

Since we are working at a fixed but large N in this paper, wc will define our length scale as 

a(b) = - 1 (33) 

It is necessary to work with a lattice coupling that shows the weak to strong coupling transition that is 
essentially free of finite N effects. We chose N = 18 and set the lattice couplings to b = 0.32,0.35,0.4. The 
fermion mass was set to fj, = 0.1 and fi — 0.05 and we set m — 4. The pseudo-fermion algorithm was used to 
compute the fermionic force. 

This is an irrelevant parameter but needs to be chosen in a specific range to realize the correct continuum 
limit. Based on previous studies QT], we set m = 4 in this paper. 

We restrict ourselves to square Wilson loops and this enables us to extract a length scale using linear inter- 
polation. The results are plotted in the right panel of Fig. [3] and compared with the two loop result for / = 
and / = 1. The length scale changes very little in the range of coupling we have studied in this paper. Let us 
assume a beta function of the form 



(3(a) = e + (a — a ) 2 



(34) 
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motivated by [5D]. Assume e > but small for our case since the two loop beta function has a zero if / is 
slightly bigger than unity. The scale as a function of the coupling is given by 

In a 2 = —= tan -1 =A (35) 

If e is small and our lattice coupling is larger than a and not close to it, the scale will change very little if we 
change the coupling. This leads us to speculate that the single site model we are simulating might be close to 
a situation where the beta function has a zero. A careful analysis of the large N corrections along with results 
at weaker coupling are needed to confirm this speculation. In addition, it will be necessary to study the chiral 
limit. The effect due to fermion masses in the right panel of Fig. [3] are small. 
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